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Recently there has been tremendous increase in the number of identified extra-solar planetary 
systems. Our understanding of their formation is tied to exoplanet internal structure models, which 
rely upon equations of state of light elements and compounds like water. Here we present shock com- 
pression data for water with unprecedented accuracy that shows water equations of state commonly 
used in planetary modeling significantly overestimate the compressibility at conditions relevant to 
planetary interiors. Furthermore, we show its behavior at these conditions, including reflectivity 
and isentropic response, is well described by a recent first-principles based equation of state. These 
findings advocate this water model be used as the standard for modeling Neptune, Uranus, and "hot 
Neptune" exoplanets, and should improve our understanding of these types of planets. 

PACS numbers: 96.15.Kc, 62.50.-p, 64.30.-t 
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The past several years have seen a virtual explosion 
in the number of extra-solar planets discovered. Two 
rapidly growing populations of exoplanets are ice giants 
referred to as "hot Neptunes" and "mini-Neptunes" ; [l| 
planets roughly the same size as or, respectively, smaller 
than Neptune and Uranus that transit their host stars at 
significantly smaller radii, resulting in higher tempera- 
tures than the ice giants in our solar system. Understand- 
ing of the composition and formation of these planets, 
and thus development of these planetary systems, relies 
on our knowledge of the equation of state (EOS) of light 
elements and compounds like water, over a wide pressure 
and temperature range. To date much of the modeling of 
ice giants has employed the ANEOS Q and Sesame H 
models for water that were developed decades ago 0, HI] ■ 
Discrepancies between these EOS models lead to signifi- 
cant differences in predicted radius evolution of Neptune- 
mass planets. Depending upon the total amount of heavy 
elements, and their distribution within the planetary in- 
terior, the resulting variation in predicted radius at a 
given age due to the water EOS can range between 5 
and 30% (f|. This is a major factor in preventing ac- 
curate determination of exoplanet internal composition 
from their observed radius. 

Recent quantum molecular dynamics (QMD) calcula- 
tions of water 0, 0| suggest an EOS that differs signifi- 
cantly from ANEOS and Sesame. Notably, when incor- 
porated into planetary models, this first-principles (FP) 
based EOS predicts a ~20% cooler core temperature for 
Neptune and Uranus The conductivity properties of 
this FP model are also noteworthy [loj], suggesting that 
water is super-ionic [ill, EH at high densities, p, and low 
temperatures, T, relevant to planets such as Uranus and 
Neptune. This predicted property plays a key role in 
dynamo models to explain the enigmatic magnetic field 



structure of these planets 13J, [14J • Another important re- 
sult is derived from the predicted phase diagram of water: 
the icy giants Uranus and Neptune perhaps contain no 
"ice" but dissociated water at a high ionic conductivity, 
even less so would close-in exoplanets. Hot and mini- 
Neptunes may even comprise water plasma with substan- 
tial electronic conduction. However, the FP EOS for wa- 
ter has not been widely accepted due to its inability to 
reproduce results from laser driven shock wave experi- 
ments in the Mbar regime [l5| . 

We present results of magnetically accelerated flyer- 
plate experiments on water performed at the Sandia Z 
machine [l6j . a pulsed power accelerator capable of pro- 
ducing extremely large current (~20 MA) and magnetic 
field densities (~10 MG) within a short circuit load. 
These data, in the range of 100-450 GPa along the Hugo- 
niot - the locus of end states achievable through com- 
pression by large amplitude shock waves - have consid- 
erably higher precision than data obtained with previ- 
ously used methods, and support the FP EOS for wa- 
ter. The high precision stems from the ability to per- 
form well-defined flyer-plate experiments on Z; the mag- 
netic pressure (>500 GPa) can propel the outer anode to 
velocities approaching 30 km/s, enabling high-precision, 
plate- impact EOS measurements in the TPa regime [TtI — 
Il9j . Furthermore, and more significantly, the present 
work obtained re-shock data of water in the range of 
200-700 GPa. These data, at high p and low T, pro- 
vide a stringent test of the isentropic response of water 
in the several Mbar regime, which is directly relevant to 
the conditions of interest for planetary modeling of Nep- 
tune, Uranus and presumably water- rich exoplan- 
ets such as the hot Neptune GJ436b 0, [H ED] . Finally, 
reflectivity on the Hugoniot was measured and compared 
to FP calculations for water I22ll. 



2 



o. 

C3 



600 
550 
500 
450 
400 



f 350 



300 
250 
200 



1 1 

Quartz drive 
plate pressure 


1 


(P \ QuartZ 

\'i' U p2)/[ \ release 






Reflected 
Quartz 
Hugoniot 


/ Quartz N 
Hugoniot 


\ 




(/••'' )\ - 

N 

1 



8 



10 12 
Particle Velocity (km/s) 



14 



FIG. 1. P — u p diagram for the water experiments for the 
case with the additional quartz drive plate. 



An aluminum flyer-plate [23[ was magnetically accel- 
erated to peak velocities of 12-27 km/s across a 3-4 mm 
vacuum gap [13] • The flyer-plate velocity was monitored 
throughout the entire trajectory using a Velocity Inter- 
ferometer System for Any Reflector (VISAR [2J),at lo- 
cations above and below an aluminum water cell [23j . A 
rear quartz window in the cell provided optical access to 
the sample. In some cases an additional quartz plate was 
placed between the aluminum drive plate and the water 
sample, enabling data to be obtained using two different 
materials as the high-pressure standard, thereby increas- 
ing confidence in the measurements. Impact with the 
cell generated a strong, multi-Mbar shock wave in the 
aluminum drive plate. This shock was then transmitted 
either directly into the water sample, or into a quartz 
plate and then into the water sample. Upon reaching 
the rear quartz window, the shock was transmitted into 
the window and reflected back into the water, which re- 
shocked the water to a higher P and p. In all cases the 
shock waves in the water and quartz were of sufficient 
amplitude that the resulting shocked material was re- 



flecting [la, LHJ |25j , enabling the shock velocities to be di- 
rectly measured using the VISAR. A total of 18 diagnos- 
tic channels were utilized for each experiment, enabling 
multiple, redundant measurements to be made, resulting 
in an overall uncertainty in the measured flyer-plate and 
shock velocities of a few tenths of a percent [231 ] . 

The shocked state of the water was determined us- 
ing the impedance matching technique and the Rankinc- 
Hugoniot (RH) jump relations [26j . a set of conditions 
derived by considering conservation of mass, momentum, 
and energy across a steady propagating shock wave. The 
shocked state of the aluminum (quartz) drive plate was 
determined from the known Hugoniot of aluminum [27j 
(quartz [191 ]) and the measured flyer-plate (quartz shock) 



velocity; this defined a point in the pressure - particle ve- 
locity (P — Up) plane, as shown in Fig.Q] When the shock 
transits into the water, a release wave propagates back to- 
ward the flyer-plate, and thus the state of the drive plate 
is constrained to lie on a release adiabat from this point 
in the P — u p plane, shown in Fig. Q] as the green line. 
The shocked state of the water is constrained to lie along 
a chord in the P — u p plane with slope given by the prod- 
uct of the measured shock velocity of water, U sw , and 
the known initial density. The intersection of these two 
curves provides P and u p , shown in Fig. [T] as (Pi,u p i); 
The RH jump relations then provide p in the shocked 
state. Uncertainties in all kinematic values were deter- 
mined through a Monte Carlo technique, which uses a 
statistical process for propagation of all random measure- 
ment errors and systematic errors in the standards p3| . 
Using this technique, the one-sigma uncertainties in P 
and p were found to be 0.5% and 1%, respectively. 

A total of 8 Hugoniot experiments were performed over 
the range of 100 to 450 GPa. Results of these experiments 
are shown as the red symbols in Fig. EJa). Also shown 
are Hugoniot data of Mitchell and Nellis [28|, Volkov 
et al. 29j, Celliers et al. [IB] ], and Podurets et al. [3(1, 
and the predicted Hugoniot response from ANEOS [2J, 
Sesame 7150 and the recent FP EOS model of French 
et al. @, [If . Note that a reanalysis of the nuclear driven 
datum of Podurets et al., using an improved aluminum 
standard for impedance matching 27J , resulted in a slight 
decrease in p. The low-P end of our data is in good agree- 
ment with the gas gun data of Mitchell and Nellis and 
the explosively driven shock data of Volkov et al. In con- 
trast, our data are significantly less compressible than 
the laser driven data of Celliers et al., which tend to sup- 
port the much more compressible ANEOS and Sesame 
Hugoniots, albeit with significantly large uncertainty and 
scatter. The vastly reduced uncertainty in p for our data, 
roughly an order of magnitude, strongly suggest that 
water is much less compressible than the ANEOS and 
Sesame models predict, and that water is instead very 
accurately described by the FP EOS of French et al. Fur- 
thermore, the reanalyzed Podurets et al. datum is also 
in very good agreement with the FP EOS. Thus, with 
the exception of the Celliers et al. data, the FP based 
model for water matches all experimental Hugoniot data 
up to 1.4 TPa. 

In all 8 of the Hugoniot experiments described above, 
the reflected shock from the rear quartz window drove the 
water from the Hugoniot state to a re-shocked state at 
higher P and p. The measured shock velocity in the wa- 
ter immediately prior to reflection from the rear quartz 
window defined the initial shocked state of the water. 
The measured shock velocity in the rear quartz window 
and the known Hugoniot of quartz provided the double- 
shocked P and Up for water, shown in Fig. [1] as (P2, u P 2). 
The velocity of the second shock in the water, U SW 2, was 
then determined by the RH jump relations using the 
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FIG. 2. (a) Water P — p Hugoniot. Models: magenta line, 
ANEOS 0|; black line, Sesame 0; dotted red line, FP 0; 
solid red line, FP 0. Data: red diamonds (squares), this work 
using aluminum (quartz) standard; open squares, Celliers et 
al. [la ]: black circles, Mitchell and Nellis [231; black triangle, 
Volkov et al. [2l|; blue (cyan) triangle, Podurets et al. [30| . 
as reported (reanalyzed [23(] ). Planetary adiabats: green line, 
Neptune (Uranus similar); blue line, GJ436b [la, H3| . 
(b) Water double-shock Hugoniot. Models: Thick (thin) red 
lines, FP first shock (re-shock) Hugoniots |2j; black lines, FP 
isentropes 0; (orange, pink, gray) lines, (FP 0, ANEOS 0, 
Sesame ,3:]) double-shock envelopes. Data: red (orange) sym- 
bols, shock (re-shock) Hugoniot data, this work. Planetary 
adiabats: as in panel (a). 



change in P, (P2 — Pi), and u p , (u P 2 — u p i). The re- 
shock p was then determined from U SW 2, the first shock 
p, and (u p 2 — u p i). Using the Monte Carlo technique, 
the one-sigma uncertainties in P and p for the re-shock 
states were found to be 0.5-1% and 1-2%, respectively. 
Although the uncertainty for the re-shock data is larger 
than that for the principal Hugoniot data (entirely due 
to the larger uncertainty in the initial state), the accu- 
racy of the present data is a significant improvement over 
previous re-shock data of Mitchell and Nellis [28| (uncer- 
tainty in p of 4-14%) and the pre-compressed Hugoniot 



data of Lee et al. [31[, (uncertainty in p of 5-10%" 



The re-shock data for water are shown in Fig. [2fb), 
where first and second shock states are correlated by like 
symbols. Also shown are several FP re-shock Hugoniots 
(thin red lines) and isentropes (thin black lines) for com- 
parison Q. These re-shock Hugoniots along with the 
known Hugoniot of quartz (i~9j were used to determine the 



double-shock envelopes - the locus of end states achiev- 
able through shock and re-shock using a quartz anvil: 
FP (orange line), ANEOS (pink line), and Sesame 
(gray line). These re-shock data further confirm the less 
compressible response of water above 100 GPa. 

Note that the FP re-shock Hugoniots (red) and isen- 
tropes (black) are nearly coincident over the p range ac- 
cessed through the re-shock experiments. This is due to 
a second order contact for the Hugoniot and isentrope 
at the initial state [26J, which is most easily seen by ex- 
panding the entropy as a function of volume in a Taylor 
series. This implies that the Hugoniot and isentrope are 
very close in P and p until, at large compression, the rise 
in T associated with the irreversible shock becomes large 
enough that thermal pressures become significant. In the 
range investigated in this study, the difference in T be- 
tween the re-shock Hugoniot states and the isentrope at 
the re-shock p, as determined by the FP EOS Q, ranged 
from 200K (out of 6800K) to 330K (out of 40000K) at 
the lowest and highest P, respectively. This makes such 
a re-shock measurement the best possible test of the isen- 
tropic response of the EOS model in this range of P and 
p. Thus the present data validates the isentropic response 
of the FP EOS in the P and p regime that is intersected 
by the water-rich models of Neptune and Uranus [a 14J , 
shown in green, and the exoplanet GJ436b 2^, 2M, 23j . 
shown in blue. 

The VISAR was also used to infer reflectivity, R, of 
water (at 532 nm) along the Hugoniot. A quadrature 
VISAR was used for all experiments, which provides 
four measures of the interference signal at 90 inter- 
vals. The signals at 180 intervals can be subtracted, 
ensuring the remaining signal only includes coherent re- 
flected laser light (incoherent light, such as self-emission 
from the hot plasma, would equally contribute to all four 
quadrature signals). Comparison of the magnitude of 
these subtracted signals before and after shock break- 
out from the water to the quartz rear window provides 
a relative measure of the shocked water R with respect 
to shocked quartz [25|. The uncertainty in R was taken 
to be the linear sum of the standard deviation of the 
inferred R from the nine independent VISAR signals ob- 
tained from each water cell and the reported uncertainty 
in R of shocked quartz [25| . 

R data along the Hugoniot are shown in Fig. [31 Also 
shown are data from Celliers et al. [l5| and the predicted 
R from FP calculations of French and Redmer [22| using 
both the Perdew, Burke, and Ernzerhof (PBE) and Heyd, 
Scuseria, and Ernzerhof (HSE) functionals for exchange 
and correlation. It was anticipated that the HSE func- 
tional, which includes the nonlocal Fock exchange, would 
prove to be more accurate in the calculation of R, as this 
functional has been shown to better reproduce the band 
gap in semiconductor materials (PBE is known to signifi- 
cantly underestimate the band gap). In comparison with 
R data of Celliers et al. 15j it would appear that the 
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FIG. 3. Reflectivity (532 nm) along the principal Hugoniot of 
water. Models: red (magenta) line, FP calculations with HSE 
(PBE) functionals Data: black line, Celliers et al. [l5| ; 

red diamonds, this work. 



HSE calculations are less accurate. However, our data 
suggest a much lower peak R, which is in significantly 
better agreement with the HSE calculations. We note 
that two recent data points (~140 and 260 GPa) from a 
group [32| at the Gekko laser in Japan also suggest lower 
R, in very good agreement with our results. These new 
results lend confidence to the FP calculations, which also 
predict a super-ionic phase of water at low T and high p 
conditions relevant to planetary interiors. Furthermore, 
these results strongly suggest that at these conditions wa- 
ter is in a plasma phase, which would imply that a T—0 
K EOS for water is not sufficient for modeling of hot and 
mini-Neptunes, and that water would be expected to mix 
in the H/He envelope rather than form an ice shell sep- 
arate from an outer H/He envelope. 

We presented data with unprecedented accuracy for 
shock compression of water to 0.7 TPa and 3.8 g/cc in a 
regime relevant to water-rich models of Uranus, Neptune 
and the exoplanet GJ436b. The experimental P, p, and 
R are in excellent agreement with density functional the- 
ory predictions, thereby validating first-principles ther- 
modynamic calculations as a sound basis for planetary 
modeling, and strongly advocating the FP EOS be the 
standard in modeling water in Neptune, Uranus, and 
"hot Neptune" exoplanets. In particular this work sup- 
ports the prediction of a ~20% cooler core temperature 
for Neptune and Uranus Q . As the calculated amount of 
H and He in the planets decreases with the stiffness of the 
water EOS, confidence in the presence of a few percent 
H and He in the deep interior of Neptune and Uranus, 
as derived from the (rather stiff) FP EOS based models 
0,0, is strengthened by this work. As H would be metal- 
lic, this might influence the generation of the magnetic 
field. Furthermore, the validation of the FP EOS in the 
regime relevant to planetary interiors all but eliminates 
one significant source of uncertainty in the predicted ra- 
dius evolution of Neptune-mass planets within assumed 
composition models. This will improve our understand- 
ing of the interior structure of these planets, and perhaps 



our understanding of these planetary systems. 
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